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Abstract. We investigate the influence of noise on a graph state generation 
scheme which exploits a mirror inverting spin chain. Within this scheme the 
spin chain is used repeatedly as an entanglement bus (EB) to create multi-partite 
entanglement. The noise model we consider comprises of each spin of this EB 
being exposed to independent local noise which degrades the capabilities of the 
EB . Here we concentrate on quantifying its performance as a single-qubit channel 
and as a mediator of a two-qubit entangling gate, since these are basic operations 
necessary for graph state generation using the EB. In particular, for the single- 
qubit case we numerically calculate the average channel fidelity and whether the 
channel becomes entanglement breaking, i.e., expunges any entanglement the 
transferred qubit may have with other external qubits. We find that neither 
local decay nor dephasing noise cause entanglement breaking. This is in contrast 
to local thermal and depolarizing noise where we determine a critical length and 
critical noise coupling, respectively, at which entanglement breaking occurs. The 
critical noise coupling for local depolarizing noise is found to exhibit a power- 
law dependence on the chain length. For two qubits we similarly compute the 
average gate fidelity and whether the ability for this gate to create entanglement 
is maintained. The concatenation of these noisy gates for the construction of a 
five qubit linear cluster state and a Greenberger-Horne-Zeilinger state indicates 
that the level of noise that can be tolerated for graph state generation is tightly 
constrained. 
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1. Introduction 

Entanglement appears to be a crucial ingredient for the potentially remarkable 
speedup of a quantum computer compared to that of a classical computer [H [2]. 
This observation is especially highlighted within the one-way quantum computing 
model [21 [4] . Here the state of a quantum many-body system, typically composed of 
spin-i or qubit subsystems, can serve a universal resource for quantum computing 
in which the computation is driven by successive von-Neumann measurements on the 
individual constituents. While the characterization of multipartite entanglement in a 
general quantum many-body state remains an open problem, initial states which can 
act as a universal resource for one-way quantum computing are within an increasingly 
well-studied class called graph states [5l [6] . 

Graph states are many-body quantum states which have an intuitive 
representation in terms of mathematical graphs. More precisely, vertices of a graph 
are assigned to the constituent qubits, each initialized in a state | = (| 0) -(- 1 1) )/ \/2, 
and edges connecting vertices represent a pattern of Ising-type interaction^ that have 
subsequently taken place between these qubits. In this way the graph describes a 
preparation procedure for this class of states, as depicted in figure [2a). Within the 
graph formalism many of the properties of graph states, such as their Schmidt measure 
and robustness to noise, can be computed efficiently despite being intractable for a 
general state [3 [71 [5] . From such studies it is known that there are graph states that 
contain the maximum amount of entanglement permissible for any given number of 
qubits. As such, graph states form a highly non-trivial class of quantum states. 

It is believed that some of the tremendous challenges faced in realizing a quantum 
computer can be lessened by using an architecture based on graph states [6]. In 
particular the underlying resource for one-way quantum computing is a special class 
of graph states, called cluster states [8^, which are represented by graphs with a regular 
lattice geometry like that shown in figure[ljb). This pattern of nearest-neighbour Ising 
interactions is a geometry which is very naturally suited to quantum lattice systems. 
Additionally, by separating the preparation of entanglement from its consumption 
within a computation the one-way model can be arranged to accommodate lossy or 
even probabilistic processes during the preparation phase. Beyond cluster states more 
general graph states are also an efficient resource for specific quantum computations [4] 
and so represent a preferred experimental route to quantum information processing 
where qubits are a precious quantity. Graph states also play a prominent role as code- 
words in quantum error correction ^ which permit the reliable storage of quantum 
information in the presence of noise. 

There arc now a diverse range of proposals for the preparation of graph states 
in realistic physical systems [6^. These include the direct use of linear optics and 
photon resolving measurements to construct graph states with photons via a non- 
deterministic protocol [To]. As a proof of principle an entirely optical creation of a 
4 qubit graph state was recently realized and used to implement a 2 qubit Grover 
search algorithm [11] . Other frameworks include using hybrid systems which combine 
optical and solid state qubits 12J. Another method is to instead engineer a many- 
body quantum system whose ground state is a graph state so that beyond engineering 
the nearest- neighbour interactions the preparation becomes a cooling problem |13j . 
The approach which we consider in this paper is based on exploiting a spin chain 

I This Ising interaction is typically taken to implement a controlled phase or c-a^ gate. 
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Figure 1. (a) An arbitrary graph state, (b) A 3D cluster state, (c) The spin- 
ladder arrangement used in the graph state generation scheme. One leg of the 
ladder is the EB spin chain with a mirror inverting Hamiltonian Hg- The other 
is a chain of decoupled spins which form a storage register R. Coupling between 
adjacent spins in EB and R is dynamically controlled to implement a rapid swap 
gate t/awap. (d) The EB spin chain with each spin exposed to an independent 
local environment £. 



with fixed engineered couplings chosen such that its dynamical evolution is mirror 
inverting [T31 [HI [IS] ■ Such spin chains have attracted much attention because of their 
ability to perform perfect state transfer and therefore act as a quantum communication 
channel [17l [TH [191 [20] . In reference [21] it was shown that mirror inverting spin chains 
are capable of implementing a specific type of multi-qubit circuit that is naturally 
suited to the generation of entanglement of the type present in graph states. For this 
reason we call this type of chain an entangling bus (EB). When the EB is used within 
a spin-ladder arrangement, as shown in figure [Ijc) where the second leg of the ladder 
is a register R of qubits, it permits the efficient generation of arbitrary graph states 
within this register. 

Experimental realizations of quantum systems inevitably possess a coupling to a 
surrounding environment composed of a large number of degrees of freedom which are 
beyond the experimenters control [22j . This coupling introduces quantum noise that 
destroys quantum coherence of the system (i.e. decoherence) . This is broadly classified 
as dissipation, when accompanied by the exchange of energy between the system and 
environment, or dephasing when there is no energy exchange. The effects of noise on a 
spin chain used as quantum channels has been investigated previously [23l [24l [23 [26] . 
Here we consider a broader set of properties including the ability of mirror inverting 
chains to both distribute and generate entanglement which are crucial for the more 
challenging use of them as EB. To do this we consider a specific, but physically 
relevant [27l [7], noise model where each spin in the EB is weakly coupled to an 
independent environment £ and the complete chain is described by a master equation, 
as illustrated in figure [IJd). 

This paper is organized as follows. In section [2| mirror inversion in spin chains is 
thoroughly described. In section[3|the graph state generation scheme exploiting mirror 
inversion is briefly reviewed. Section [4| outlines the methods we apply to characterize 
the performance of the EB as a single-qubit channel and as a two-qubit gate in the 
presence of noise. In section [5] the class of local noise that is considered in this 
work is introduced. The influence of these local noise models are then systematically 
analyzed in section [7| for both the single-qubit channel and two-qubit gate scenario. 
For the readers convenience the detailed results of section [7| are summarized in its first 
subsection. We then examine the implications of these results for the generation of a 
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five qubit linear cluster state and Greenberger-Horne-Zeilinger (GHZ) state in section 
[5] before concluding in section [31 

2. Mirror-inverting spin chains 

Our starting point is a spin-i chain composed of N spins which is governed by an XX 
Hamiltonian of the form (taking h = \) 

N-l N 

with spatially dependent spin couplings tj and local fields hj . We denote the cr^ basis 
states of the chain as \qi, . . . , qN) with qj G {0, 1} representing "f and J, respectively. 
Since [Hs,¥l] = 0, where N = ^EjloCl ~ ^|) B ^hen Hs is block-diagonal with 
respect to subspaces Tin spanned by states \qi, . . . , qN) with J2j Ij — ifi- The spin 
chain Hamiltonian Hs can be mapped to a ID spinless fermionic lattice model using 
the Jordan- Wigner transformation (JWT) [28] giving 

JV-l N 

hf= -jYI + sVi^j) + E '^j^^j' 

and subsequently diagonalized into an explicitly free-fermion bi-linear form 

N 

Hd = E 

k=l 

with spectrum e^. Here both cj(cj) and a\,[ak) are fermionic creation (annihilation) 
operators, obeying the usual anticommutation relation, associated to lattice site j 
and the energy eigenstate k respectively. Under this mapping N — cjcj and the 
subspaces Tin it defines are identified with the fermion number. The fermion vacuum 
is then | vac) = | 0, . . . , 0) with energy E'vac = and spin states | gi, . . . , q^) become n 
fermion Fock states | gi, . . . , qjyi) ^ (c|)^^ . . . (cjy)'^" | vac) with the operator ordering 
following the lattice numbering. We denote the blocks of Hp acting on subspaces 
Tin as H^'p^ and since Hp is a non-interacting Hamiltonian its properties are entirely 
defined by its single-particle Hamiltonian H'^p^ . 

To be mirror inverting all localized states | j) — Cj \ vac) in Tii are required to 

evolve after a given fixed time r under Hp^ into the localized state | j) (up to a 
phase) where j = iV — j -I- 1 is the mirror location in the lattice. While this places 
constraints on the couplings tj and fields hj there are still an infinite number of 
permissible choices [T7| [H]. In this work we exclusively consider the simplest and 
fastest mirror inverting couplings [19] where tj = j{N — j) and hj = h. With this 

choice h'^P takes the form = —JSx + h where Sx is the x-axis angular momentum 
operator for a spin-5 pseudo-particle where S = ^{N — 1). Localized n = 1 states are 
then identified with Sz eigenstates {\S,l)^} of the pseudo-spin through the ordering 
1 1) = |5, — 5)^, . . ., I A^) = I 5,5)^. If we now consider the time evolution in Hi 
for a time t — n/J we see that C/'^-' — exp(— iijj^,"'^V) — exp(— i/itt/ J) exp(i7r52;) is a 
rotation of the pseudo-spin by tt about its z-axis and is therefore equivalent to the 

§ The operator N counts the number of spins which are J,. 
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C{N) 

Figure 2. The quantum circuit C{N) composed of c-ct^ gates between all 
distinct pairs of qubits obtained by evolving the mirror-inverting spin chain with 
Hamiltonian for a time r. 



mirror inversion of a single fermion in the lattice |14[ 115] . Interestingly, we note that 
similar effects also take place in a chain of coupled harmonic oscillators as discussed 
in reference [29]. 

Moving our consideration back to the full state space of the lattice it follows 
that the localized modes cj are related to the energy eigenmodes a\. via irreducible 
representations djk (f ) of a 7r/2 rotation about the y-axis of the pseudo-spin [30] as 
a\ = djk (-1) Cj. The angular momentum couplings also result in the spectrum of 
Hu being linear as = J{k — S — l) + h over the range efc G [—S + h,S + h], and so in 
order to ensure that the state | vac) is the non-degenerate ground state of the system, 
for all J, we require h > SJ. We can define the many-body gap between the vacuum 
ground state and the first excited state as A giving h = SJ + A and for A > the 
first excited state is always in the Hi subspace. Note also that with this definition in 
the limit J ^ we have that A is the local gap for each decoupled spin. Additionally, 
we can choose A/ J as an even number such that mirror inversion proceeds with no 
phase modulo 2tt. 

The mirror inverting dynamics in Tii is equivalent to the transformation UcjU^ = 
ct on the localized modes with U = exp{—iHpT). Applying this evolution to an 
arbitrary n fermion Fock state, and performing the inverse JWT, mirror inversion 
results in 

where !]„ — in(n— 1) is the number of anti-commutations of the operators cj required 
to reestablish the correct ordering. The simplest utilization of mirror inversion is state 
transfer where we restrict our consideration to the subspace Ho Hi spanned by the 
spin-polarized state | vac) and the single spin-flip states | j). We then encode an input 
qubit as a superposition \ ip) = fo | vac) -|- i^i 1 1) using the first spin in the chain and 
under purely coherent evolution this state is transferred perfectly to the last spin as 
\tp) = fo I vac) + vi I N) [HI [T5] . The same conclusion follows trivially for a mixed 
input state. 

A more general use of mirror inversion follows from noting that the phase 
TrEn in equation ([2]) is non-linear in n and only appears between subspaces with 
different fermion number for n > 2. Thus for input states of the chain which 
involve superpositions spanning several multi-particle subspaces these phases will 
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create entanglement in the mirror- inverted output state [31] [T7j . More precisely, the 
evolution U of the chain for a time r is equivalent to a quantum circuit C{N) composed 
of c-cr^ gates between all distinct pairs of N qubits followed by the inversion operator 
M, as shown in figure^ This circuit has the useful property that if any N — q spins in 
the chain are in the state | 0), then this circuit reduces to C{q) between the remaining 
q qubits, independent of their locations, followed by the full inversion M of the chain. 

3. Graph state generation with an engineered spin ladder 

Here we briefly review the scheme given in |21) where the general multi-qubit circuit 
C{N) implemented by a mirror inverting chain is exploited to construct graph states. 
This is achieved by considering a spin-ladder with a comfe- like arrangement of couplings 
as depicted in figure [Uc). One chain of the ladder possess fixed mirror inverting 
couplings and forms the EB, while the other chain is composed of decoupled spins 
forming the register R. We assume that spins in the register can be individually 
manipulated and measured. Dynamical control of the spin couplings is restricted 
to those between adjacent spins in the EB and R where we require the ability to 
rapidly implement a swap gate. In this way entanglement generation is achieved by 
repeatedly swapping qubits between R and EB and thereby using the quantum circuit 
C{N). 

The entire spin ladder is taken to be initialized in a spin polarized state. The 
scheme begins by choosing a set of register spins G that will be the graph qubits, and 
transforming all of them to | . For any subset Q C G of graph qubits which are 
transferred into the EB and evolved for a time t the resulting circuit C(|(5|) will apply 
c-a^ gates between all of the corresponding graph vertices. In the case where two 
graph qubits in the set Q do not already possess an edge between them this process 
will establish one, otherwise it will remove the edge. By proceeding iteratively we can 
induce any pattern of edges between the graph qubits G. Starting with g — I, we 

(i) transfer the g-th graph qubit from G, and all graph qubits gc > g which are 
required to connect to g, as specified by the graphs adjacency matrix F, into the 
EB; 

(ii) allow the EB to evolve for a time r and create a complete set of connections 
between all these previously unconnected vertices; 

(iii) then transfer qubit g back to the register while leaving the qubits gc to evolve for 
one cycle longer in the EB, subsequently removing all the connections between 
them; 

(iv) finally the qubits gc are transferred back to the register and step (i) is repeated 
with g 1-^ g + I. 

Thus, any graph with n vertices can be generated in at most 0{2n) uses of the EB 
in contrast to 0{n^) steps if the EB was used to implement single c-a^ gates only. 
Although the EB has a linear topology, by using this method any two qubits in the 
register can be entangled thereby allowing for arbitrary topologies of the graph state. 
To avoid overlap between EB and register graph qubits after inversion one may choose 
|G| < with locations in the first half of the register. 
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4. Characterizing a noisy spin-chain 

The main aim of this work is to characterize the effect of noise on the performance 
of the EB and determine its imphcations for using the EB within the graph state 
generation scheme. For simplicity we determine the performance of the EB at 
implementing its two most basic operations, namely acting as a single-qubit quantum 
channel and as a mediator of a two-qubit c-cr^ gate. These represent the minimal 
operations required for the EB to be used for graph state generation. For this reason 
we consider the effect of noise only on the EB spin chain and not the register R. 
Additionally we focus on the scenario in which the input and output qubits are the 
end spins of the EB. Before describing any specifics about the noise we first outline 
some general theoretical tools which provide insightful measures of performance. 

^■l. Average fidelity 

Suppose we have a system which, when no noise is present, performs a particular 
unitary operation U . With the inclusion of noise the action of the system is instead 
described by a superoperator A. How close the noisy operation remains to U for a 
particular initial pure state | V') G can be quantified by the fidelity [31] 

F{i,) = {,k\U^ {^\}U\^) . (3) 

The overall performance of the noisy system at implementing U can then be measured 
by the average of this fidelity over all possible initial pure states 

{F)= [ F{^)d^lj, 

where integration is over the unit sphere S'^'d-i in (C^ and dtp is the normalized 
measure on the sphere, also known as a Haar measure. For the case of a 
single qubit this is equivalent to integration over the Bloch sphere as J^a dip = 
J^^d(t) Jq d9 sin(6'). Now given a Kraus decomposition of the superoperator A 

as 

d= 

m— 1 

where Am are Kraus operators there is a compact formula for (F) in any dimension 
d. Firstly, we form a new superoperator £ with Kraus operators Em = AmU\ such 
that £{U pW} = A{p}, which describes exclusively the effect of noise. It can then be 
shown [31 [33] that 

(^> = ^(E|tr(i.„OP + dj. (4) 

We exploit this formula to determine the single qubit channel (or 1 operation) fidelity, 
and the gate fidelity for the effective c-tr^ operation between two qubits achieved with 
a noisy EB. 

4-2. Entanglement breaking and generation 

While the average fidelity provides a quantitative measure of a noisy operation, a more 
qualitative way of characterising the severity of the noise is to determine whether 
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Figure 3. A schematic representation of the Jamiolkowski isomorphism used 
to characterize a superoperator (a) A^^l acting on a single subsystem and (b) 
A[21 

acting on a pair of subsystems in the maximally entangled state | 't) (see 
[Appendix C| for details), as quantum states and (A respectively. 



the corresponding superoperator A^-'^l, which acts on one subsystem, preserves any 
entanglement that the subsystem has with other external systems. Quite generally 
if AI^I acts on the subsystem 6, with Hilbert space Tij, = C^', it is described as 
entanglement breaking if the final state — la ® ^b^{p^ab\ separable for 
every (possibly entangled) initial state p™^^ of the composite system of h and another 
subsystem a with Hilbert space Tia — C"^" . Becoming entanglement breaking therefore 
signifies that the channel can no longer be used to distribute entanglement. 

Remarkably, for a singlc-qubit (db = 2) the PPT criterion [331 [Mj (see 
[Appendix B[ ) in combination with the Jamiolkowski isomorphism [36| (sec |Appcndix C"| 
and figure [3]) give a straightforward condition for A^^l to be entanglement breaking. 
Firstly, it is sufRcient to compute the state from the Jamiolkowski isomorphism 
(see figure [3l[a)), where da = <db = 2, since this contains all the properties of A^^l. It 
then follows that A^^l is entanglement breaking (for any da) if and only if the state 
is separable since this implies that A^^l has a Kraus representation composed entirely 
of projectors. Finally, since p^ is a two qubit state its separability follows directly 
from the PPT criterion. The entanglement breaking characteristics of the EE when 
acting as a single-qubit channel are of importance since the graph state generation 
scheme involves its successive use. We therefore have a minimum requirement that 
for the EB to be useful it must, at the very least, preserve any entanglement that an 
input qubit may have with other external qubits, such as those in the register, when 
acting purely as a quantum channel. This then provides an essential, albeit optimistic, 
bound to its tolerance for noise. 

When the evolution of two subsystems is described by a superoperator A^^l 
it is clearly of interest to determine when this evolution is capable of generating 
entanglement between these subsystems [[j. Specifically, for pg"* = ^ic'i/^bc} '^^Y 
ask when is p™' always separable for all possible separable initial state p™^l This 
implies that the superoperator A^^l never generates entanglement. The question can 
be answered by again appealing to the Jamiolkowski isomorphism via the state 
associated to Al^l. It follows that A^^l is of product form A|j^' — Aj^^' (g) A^^' and 
incapable of generating entanglement if its corresponding state is separable with 
respect to the bipartition of the system as {ah){cd) as in figure [3jb). Thus the property 
of entanglement generation can also be phrased as a state separability problem. For 
two qubits (da = d;, = 2) the mixed state describes four qubits. In this case the 



{{ When this results in an entangled mixed state it can then, in principle, be distilled. 
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PPT criterion only provides a necessary condition for the {ah){cd) separability of this 
state. Thus the PPT criterion can only determine a point at which we can no longer 
be certain whether A^^l can generate entanglement. Nonetheless this point provides a 
quantitative cut-off which should be avoided if the noisy entangling operation is to be 
of practical use. 



5. Noise models 



We consider noise which is described by a quantum master equation of Lindblad form 
^^p{t)= -x[Hs,p{t)]+C{p{t)}, (5) 

where £{•} is the Lindbladian describing the incoherent contribution to the evolution 
of the density matrix p{t). The microscopic derivation of such a master equation relies 
on the Born-Markov approximation and is typically found to be accurate for systems 
with a weak coupling to a much larger environment [22] . 

We consider a subclass of this noise model where each spin experiences an 
independent local environment so the Lindbladian decomposes as a sum £{•} = 
We make one further restriction and consider the local Lindbladian Cj{-} 
to be of a physically well motivated form commonly encountered in quantum optical 
problems after applying the rotating wave approximation [27) [J • Specifically, 

I [2a+p{t)aJ - aja+p{t) - p{t)a- a+] 

[a|p(<)a|- pW], (6) 



2 

-I 

where a, (3 and 7 are the rates for jumps 1 1) | i), | i) IT): and pure dephasing, 
respectively. To give an overview of the physics contained in this model let us consider 
the situation where J = in iJg, defined in equation ([1]), so each spin decouples with 
a local Hamiltonian of the form Hj — -^(1 — o"!). If we move to the interaction picture 
of Hj for each spin0 it is readily found that the evolution of the j-th spin is described 
by p(t) = e^^^{p] with Cj remaining in terms of the untransformed operators in 
equation ^ due to phase cancellation. We now write the initial state p as 

p= \ + {ax)dx + {a-y)cFy + {az)crz\ , 
from which the general solution is found to be [27] 



As is well known this solution shows exponential convergence with rate a -1-/3 of {az(t)) 
to its stationary {t 00) value of (cr^)^ — and the exponential decay, with rate 

^{a + (3 + 27), of the coherences (a±(t)) to their stationary value (cr±)s = 0. The 
general solution to this noise model can be expressed in a Kraus form [31 with Kraus 

^ We shall denote the interaction picture of a Hamiltonian by a tilde as O ■■ 
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- ' T2 j ' ^2 - V n 

TsOA^ rrr/ 



operators 



where P| = ' ~ stationary spin populations and T j are functions of 

a, (3 and 7 which we give exphcitly in [Appendix A[ We also show in [Appendix A| that 
this noise model reduces to a number of well known and simpler models in specific 
limits. In particular, if we parameterize the rates as (taking Boltzmann's constant 

kB = 1) 

for an arbitrary 7, we obtain a total decay rate a(T) + f3{T} ~ k that is independent 
of T. In this case the stationary density matrix ~ limj^oo p(t) for the spin is 
equivalent to a thermal state of temperature T. Consequently in this regime the 
master equation describes finite-temperature noise caused by the coupling to a generic 
thermal reservoir at temperature T local to each spin. At T = and 7 = wc have 
a = and (3 — k which describes decay noise. For T — > cxd we have a = /3 = and 
after setting 7 = ^k, so the populations and coherences decay at the same rate, this 
results in depolarizing noise. Taking a = P = and 7 = k we obtain pure dephasing 
noise. We shall consider each of these limiting cases as local noise in the EB. 

If we now consider this class of noise in the context of a single-qubit channel 
we can determine the properties which were outlined in section |4l Indeed using the 
Kraus operators in equation ((Tj) the average fidelity of the channel can be computed 
via equation ^ and is given in full in [Appendix A[ The channel can be shown to 
become entanglement breaking if and only if the following condition is satisfied [7] 

2P|Px e^''*{cosh{[a + (3]t) - 1} > 1. (9) 

In particular this result immediately indicates that, regardless of 7, whenever a = or 
/3 = the channel is never entanglement breaking for finite coupling k and times t since 
either P| = or P| = 0, respectively. This result similarly holds when both a = 
and /? = 0, giving a pure dephasing channel 0, since we have cosh([Q! + — 1 = 0. 
In contrast a finite temperature channel (for any 7) can always become entanglement 
breaking for a finite k and t. For 7 = this entanglement breaking occurs for a 
coupling 

■(1 + 6-^/^)2 



J , 

> — cosh 

TT 



taking t = t, and this saturates at Kc ~ 0.56 J for T 00. The presence of dephasing 
reduces this threshold. An important special case is where 7 = which in the 
T ^ 00 limit gives a depolarizing channel with a threshold Kc> log(3) ~ 0.35 J. 



6. Numerical method 



The numerical calculations we perform in this work is restricted to the class of ID 
quantum lattice systems described by a master equation which include a Hamiltonian 

Note in the case of pure dephasing the definitions of P| and in terms of a and /3 are meaningless 
and the stationary populations follow from the arbitrary initial state. 
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(a) exp(;&r) 
(b) 




Figure 4. The setup used in numerical calculations to determine the effective 
two qubit superoperator API of a noisy mirror-inverting chain, (a) Following the 
Jamiolkowski isomorphism the initial state p is a spin-polarized chain | ■ ■ ■ 0) 
aside from the end spins which are in maximally entangled states | ^) with 
corresponding ancillary spins. The spins 1, ■ • ■ , A'^ in the chain are then evolved 
for a time r while being exposed to noise. The total dynamical evolution of the 
chain is then described by the superoperator exp(i3r) which is the formal solution 
to equation JSJ. (b) The state corresponding to A PI is then extracted from 
the overall final state exp{BT){p} by tracing out all but the end spin pairs. 



and a Lindbladian that are both composed of terms involving at most nearest- 
neighboring sites. It can be seen that both Hs and £ introduced in section [5] 
satisfy this constraint. The real time evolution for this class of master equation can 
be computed efficiently and to near-exact precision for systems composed of many 
sites using the mixed-state version of the Time Evolving Block Decimation (TEBD) 
algorithm [371 121 13H] ■ We refer the reader to the literature for a detailed description 
of this method and note here only that for the calculations presented we found that 
a truncation parameter [38j up to x = 20 was sufficient. In figure H] we depict the 
type of numerical calculation we have performed with this algorithm. These are based 
on the Jamiolkowski isomorphism which for two qubits requires the spin chain to be 
initialized in a pure state | VE") (8) | • • • 0) (8) | 5') with the two end spins being in a 
maximally entangled state | = (| 00) + | ll))/-\/2 with corresponding ancillary spins 
shown in figure HJa). We then use TEBD to time evolve the spins 1, • • • , iV in the 
chain in the presence of noise and finally compute the reduced density matrix of 
the two ancillary spins and the two end spins of the chain as depicted in figure IH^b) . 
The state then completely characterizes the accumulative noisy operation of the 
chain A^^l for two qubits. 

7. Results 

7.1. Summary of results 

Having introduced all the necessary concepts we now investigate the influence of 
local decay, dephasing, thermal and depolarizing noise on the performance of the 
EB spin chain. We distinguish between two scenarios, namely where only one qubit is 
transferred into the spin chain so it acts as a quantum channel, and where two qubits 
are swapped into the chain such that the mirror inversion performs a c-cr^ gate. We 
summarise the main results here and refer the reader to the proceeding subsections 
for more details. 
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When used as a single-qubit channel we find that neither local decay or dephasing 
noise become entanglement breaking. For local decay noise we find that the 
accumulative noise A 1^1 of the chain is identical to the local decay noise on any single 
spin and is therefore entirely independent on the chain length N . We show this useful 
property is a consequence of commuting coherent and dissipative contributions to the 
dynamics which is unique to local decay noise. The EB chain is found to be the most 
robust to local decay noise and is able to maintain {F) > 0.99 for state transfer when 
k/J< 9.7 X 10"'^. The case of local dephasing noise is shown to be well modelled 
by a quantum channel subject to the same local dephasing noise along with a length 
dependent decay noise. For spin chains of lengths up to = 50 the average fidelity 
remains above 99% as long as the dephasing noise fulfils k/J < 5 x 10^'^. 

In contrast to these types of noise we find that local depolarizing and thermal noise 
become entanglement breaking for certain parameter regimes, for which analytical 
estimates are given. We further find that the length independence observed for the 
T = local decay noise persists as a very weak length dependence for significant 
non-zero temperatures T < 0.2 A. For the local depolarizing noise we find that the 
critical coupling at which entanglement breaking occurs for a given chain length N is 
described by a power law Kc/J « iV~^ with x = 0.68. This behaviour appears to be a 
consequence of the competition between the speed and spreading of a spin-packet in 
the chain. As expected the chain is most severely affected by local depolarizing noise 
with (F) > 0.99 only for k/J < 3 x 10"'' and lengths up to iV = 50, which is more 
than an order of magnitude smaller than for local decay or dephasing noise. 

For the two-qubit case we find that the average gate fidelity with local decay 
noise does not depend on the length TV and remains above 99% for couplings below 
k/J < 4 X 10"'^. For thermal noise we find that (F) only depends very weakly on N as 
long as T < 0.2A, as was the case a single-qubit. However, in contrast to the single- 
qubit case, this length independence in (F) does not extend to the accumulative noise 
AI^I superoperator itself. Both local dephasing and depolarizing noise have average 
gate fidelities which depend on N and for up to = 12 spins the coupling is restricted 
to k/J < 2.5 X 10-3 and k/J < 4 x lO"'^ in order for (F) > 0.99. This again 
indicates that local depolarizing noise has the most severe influence and explains why 
its accumulative noise is well approximated by product noise. 

7.2. Single-qubit channel 

In this section we consider the EB as a single-qubit channel and systematically 
compute the average channel fidelity (F) and the minimum eigenvalue emin of the 
partial transposition of the mixed state that is isomorphic to accumulative noise of 
the chain A^^l. From the PPT criterion this noise is entanglement breaking whenever 
Cmin > 0. We also use the behaviour of (F) over a wide parameter range to fit the 
noise A'^l of the chain to the specific class of single-spin noise introduced in section [5] 
and find that such fits are possible to very good accuracy. 

7.2.1. Decay noise (low-T limit) - As was noted in reference |24j . we find from 
our numerics that (F) displays no dependence on the length of the chain. Here we 
show that this unexpected feature is in fact a consequence of a much stronger result; 
specifically, the superoperator A^^l itself, which characterises the accumulative noise 
of the chain, is independent of N. This result implies that A^^l for any N is equivalent 
to AI^I for a chain with = 1. Since a = 1 chain is a single-qubit decay channel. 
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this allows us to conclude that state-transfer in a mirror-inverting chain with local 
decay noise is never entanglement breaking. Additionally, the coupling at which the 
fidelity drops to (F) < 0.99 is found to be K//J = 9.7 x 10~^, independent of N. 

We now explain the independence of chain length observed. To begin we take a 
spin chain composed of N spins and a general Hamiltonian Hg satisfying [Hs , N] = 
so Hs is block diagonal with blocks -ffi""* in each subspace Tin Q- Then we restrict our 
considerations to initial states of the chain p(0) whose support is entirely contained in 
the subspace TLq © Tii . Evolution due to Hg and local decay noise has the convenient 
feature that the support of the state pit) at any time will also remain entirely within 
Tio ffi Til • Consequently, we may project the full master equation of the chain into this 
subspace yielding 

= - ® Hi'\p{t)] + I (2Fi(t)Po - Pip(i) - p{t)T,) . 

Here Pi{t) — tr(Pip(i)Pi) is the probability of being in the Tii subspace, with Pq and 
Pi being the projectors onto the subspaces TLq and Tii, respectively. If this projected 
master equation is expressed as p{t) — H{p(t)} -I- L{p(t)}, where H and L are the 
coherent and dissipative superoperators, it follows that [H, L] = since L is composed 
entirely of projectors onto the same subspaces over which Hg is block-diagonal. The 
crucial effect of this commutivity is that 

p{t) = e"*+"{p(0)} = e''*e"{p(0)} = e"e"*{p(0)}. 

Hence the coherent and dissipative contributions to the evolution are independent and 
can be applied separately. 

In the special case where Hg — Hs is a mirror-inverting Hamiltonian this property 
manifests itself directly in the accumulative noise A^-'^l. Using the chain as a channel 
involves initializing it in a spin-polarized state 

-o(0) = I T> (T li ® • • • ® ® • • • ® I T) (T Ij ® ■ • • ® I T) (T U . 

aside from the spin j which is in input state <7. If we first apply the coherent evolution 
for a time r then, as outlined in section [21 the state becomes 

Pi = e^'^ipio)} - I T> (T li ® • • • ® I T> (T I, ® • • • ® 9 ® • • ■ ® I T> (T U > 

where the state <; is transferred to the mirror spin j. Since the whole chain is spin- 
polarized, aside from at spin j, the action of the superoperator L on such a state is 
completely equivalent to single-spin decay noise at that spin alone. Thus, the final 
output state ip of spin j after tracing out all other spins (denoted as c) 

= tr,(e^^{pi}), 

is identical, irrespective of N [|, to the output state for a chain with iV = 1 where the 
input state ^ is simply exposed to single spin decay noise for a time r. 

7.2.2. Dephasing noise - We find that the behavior of Cmin and {F), displayed in 
figure [5ja)-(b), for the accumulative noise of the spin chain with local dephasing 
rapidly converges with the length of the chain. In particular figure [Ifa) shows that 
the chain does not become entanglement breaking for the wide range of chain lengths 
N and couplings k/J investigated. In figure[5jc) we plot the coupling i^f/J at which 
the fidelity drops below {F) < 0.99. This plot indicates that the coupling must not 




Figure 5. For local dephasing noise - (a) The minimum eigenvalue tmin of the 
partial transposition of and (b) the average fidelity (F), both as a function of 
the chain length A'^ and coupling strength ft/ J. (c) The coupling Kjj J at which 
the fidelity shown in (b) drops below (F) < 0.99 against the chain length N . (d) 
The fit parameter as a function of N . 



exceed k/J ~ b y. 10^^ for chain lengths of order = 50 for useful fidelities to be 
achieved. 

One might expect that the average fidelity (F) can be reproduced by assuming 
that a single qubit is sent through a purely dephasing channel with a coupling k 
dependent on N. However, our numerical calculations show this not to be the 
case. Instead, we find that {F) is fitted extremely well by assuming that the overall 
noise A'^] is simultaneously decay and dephasing. This model noise is also never 
entanglement breaking for finite n/ J . Using the general expression for the noise model 
in [Appendix A| we fitted {F) for each N to the model noise fidelity with 7 = k and 
(5 — ^K, where C is the only fit parameter. The parameter C, obtained as a function 
of N is plotted in figure [nfd). It shows that as the chain length increases the decay 
rate (3 increases and becomes of the same order as the dephasing rate 7. Intuitively 
this type of model might be expected to describe the accumulative noise of the chain. 
Like local decay noise in the section [7. 2. II an initial state with support in 7io ©T^i will 
remain so at all times. However, unlike local decay noise the coherent and dissipative 
contributions to the projected master equation do not commute and so perfect mirror 
inversion is not obtained for any k/ J > 0. As a result the input state on the first spin 
is never perfectly refocussed to the A^-th spin causing a 'leakage' of the [ population 
over other spins in the chain. Since the A'^-th spin is the output qubit at time r this 
effect appears as decay noise. 

* The Hamiltonian Hs defined in section [2] is one such example, 
tl On the proviso that the inversion time t is kept constant with N . 






Figure 6. For local thermal reservoirs - (a) The minimum eigenvalue Emin of the 
partial transposition of and (b) the average fidelity {F), both as a function of 
the chain length TV and temperature T for noise caused by a coupling strength 
k/J = 0.02. (c) The temperature 7//A at which the fidelity shown in (b) drops 
below 99% of its value at T = against the chain length A'^. (d) The fit parameter 
as a function of N . 



To gauge how accurate the assumed noise model was compared to the actual noise 
superoperator A^^l we computed the fidelity [31] 

F^{p\pt) ^ tr [47p^i^t7p^) (11) 

between the states and isomorphic to A^^l and the model noise, respectively (see 
[Appendix C| for details). We find that over all parameters considered the infidelity 
1 — i^A < 3.2 X 10~^ which indicates that the model is capturing the accumulative 
noise of the chain to good approximation. 

7. 2. 3. Thermal noise (finite-T ) - For local thermal noise we restrict our consideration 
to a suitably weak coupling k/J = 0.02 so that the corresponding average fidelity at 
T = is {F) = 0.98 and still sizable. In figure [SJa) emin is plotted and demonstrates 
that up to temperatures T = A and chain lengths iV = 50 the accumulative noise 
of the chain is not entanglement breaking. The behaviour of both emin and (F) in 
figure [Bl^b) indicates that their inscnsitivity to the chain length N, seen earlier for the 
T = decay noise, persists for temperatures T < 0.2A. This is further confirmed by 
figure[SJc) where the temperature Tf/A at which the fidelity drops to 99% of its value 
at T = is above T — 0.2A for chain lengths up to = 50. 

In order to reproduce the fidelity surface of figure [6]^b) we fitted a noise model in 
which a{T) and /3(T) remain unchanged from those in equation ([5]) but now include 
a non-zero, temperature-dependent dephasing rate j(T) = (^a{T). As a result the 
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effective noise of the chain is still described by a coupling to thermal reservoir of 
temperature T. The corresponding fit parameter C, for each N is shown in figure [Hid) 
and is seen to be linear and very nearly C,{N) = N . The accuracy of this model noise 
compared to the numerically determined noise was found to be extremely good with 
-Fa computed via equation pT|) giving 1 — Fa<3.6x10~^. 

In section[5]we found from equation (|10p that for single-qubit thermal noise, with 
7 = and acting for a time t, a coupling k/J > 0.56 was required for the channel 
to become entanglement breaking at T — > oo. Consequently, a single-qubit channel, 
with the weak coupling k/ J — 0.02 chosen, never becomes entanglement breaking at 
any temperature. Using our model for the accumulative noise of the chain in which 
a non-zero ^{T) « Na{T) emerges we have determined an approximate analytical 
expression for the critical length Nc at which entanglement breaking will occur for a 
given temperature T and local coupling k/J as 



For any k > this function monotonically increases with decreasing T from a finite 
asymptotic value at T — > cxd and diverges at T = 0. For the weak coupling used in 
this section the critical length at T ^ oo is A^c = HI spins. 

7.2.4- Depolarizing noise (high-T limit) - For local depolarizing noise we find that 
the accumulative noise of the chain becomes entanglement breaking at a threshold 
coupling Kc/ J that reduces with the chain length N as shown in figure [7ja). The 
fidelity {F) shown in figure [Tf^b) decreases rapidly with k/J for N > 5. The coupling 
Kf/J at which the fidelity drops to (F) — 0.99 is plotted in figure [Tjc) and indicates 
that the condition k/J < 3 x 10^^ has to be fulfilled in order to achieve reasonable 
fidelities for chain lengths up to iV = 50. Figure [7]^b) also shows that the decay of (F) 
with k/ J changes from an exponential behaviour for small iV to a double-exponential 
behaviour for large N. This indicates that the accumulative noise of the chain does 
not remain purely depolarizing, however, we do find that it is still well approximated 
by the class of noise introduced in section [51 Using the noise model with a = /3 = CiK 
and J = and fitting ^1 and ^2, the fidelity curves can be accurately reproduced for 
all parameters considered. By restricting a = f3 this model is still thermal noise in the 
limit T 00, but importantly we allow the total decay rate to increase from k and 
also independently allow the dephasing rate to increase from ^k, both as a function 
of N. The fitting parameters plotted in figure \T[d) show that Ci increases from its 
initial value of ^ to a little over unity, whereas (2 displays a linear increase with N 
becoming nearly 30 times larger than Ci for iV = 50. To establish the validity of this 
model we compare it to the actual noise computed numerically by calculating F\ via 
equation (11). We find that 1 - Fa < 1.5 x lO^^. 

The dependence of the critical couphng Kc/J with N is plotted in figure [8] and 
appears to be described well by a power-law Kc/ J ~ with x = 0.68. Using the 
fitted noise model the critical coupling can be obtained by solving a special case of 
equation ^ of the form 



as a function of TV using the functions Ci(-^) ^^'^ C2{N) plotted in figure [Tljd) . The 
result of this is also shown in figure [5] and is consistent with a power law with exponent 





X = 0.72. 
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Figure 7. For local depolarizing noise - (a) The minimum eigenvalue emin of the 
partial transposition of and (b) the average fidelity (F), both as a function 
of the chain length A'^ and coupling strength k/J. The critical coupling Kc/J for 
each chain length A'^ at which the accumulative noise is entanglement breaking is 
given by the intersection with the emin = plane shown in (a), (c) The coupling 
K//J at which the fidelity shown in (b) drops below (_F) < 0.99 against the chain 
length N . (d) The fit parameters Qi (left axis and 'o') and ^2 (right axis and 'o') 
as a function of A''. 



10" 




Figure 8. The critical coupling Kc/J at which the accumulative noise of the 
chain becomes entanglement breaking as a function of the chain length A'^ on a 
log-log scale. The numerical data is plotted with 'o' and the fitted power law 
Kc/J ~ is the solid line with an exponent x = 0.68. The solution for the 

critical coupling of the noise model are plotted with 'o' and the dotted line is the 
power law fit with an exponent x = 0.72. 
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Figure 9. (a) A simple model for state transfer in a mirror inverting chain with 
local depolarizing noise. Each spin on the chain is considered to be a depolarizing 
channel T>j{p} = pjp + ^(1 — Pj)l with pj = exp(— k/jSjt). The fraction of the 
total time t spent in each channel fj is approximately proportional to the inverse 
of the average spin-coupling for the j-th spin, (b) The spin packets probability 
distribution \cj\^ in the \j) basis for a selection of times. These are readily 
computed from the x-axis rotation of z angular momentum states from which the 
spread Sj can then be extracted, (c) The fraction of time fj (left axis and solid 
line) and spread Sj (right axis and dashed line) as a function of the spin j. 



To gain a qualitative understanding of the origin of this power-law scaling of 
Kc/J with N we consider a simple model of this noise scenario. Specifically we replace 
each spin of the chain by a depolarizing channel which preserves the input state with 
probability pj and where state transfer corresponds to the concatenation of these 
channels shown in figure [H^a). The accumulative noise of this sequence of single qubit 
channels is then also a depolarizing channel with probability p = Y[j Pj ■ In figure (Hlb) 
the progression of a spin packet in the chain resulting from a spin-flip excitation at 
the first spin is shown for a sequence of times for = 50. With this in mind we take 
the probabilities for each channel as pj ~ exp{—KfjSjT) where fj is the fraction of the 
inversion time r the centre of the spin packet spends at spin j, and Sj is the number 
of sites the spin packet spreads across when it is in the region of spin j. Both these 
quantities can be readily derived from the properties of angular momentum and are 
plotted in figure [9l^c). From this we see that the spin packet is narrow and slow at the 
edges while being wide and fast at the centre. The critical coupling for this simple 
model can be extracted from equation ^ as 

Kc _ log(3) 

We find that the competition between the spreading and speed of a spin packet across 
the chain as a function of its length N naturally gives rise to a power-law dependence 
for Kc/J- 
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7.5. Effective two qubit gate 

In this section we consider the EB as a mediator of an effective two-qubit c-tr^ gate. 
We systematically compute the average gate fidelity {F) for this operation and the 
minimum eigenvalue £min of the partial transposition, with respect to the bipartition 
(12)(34), of the 4 qubit state g isomorphic to the accumulative noise < 

then the noisy operation of the EB is still capable of entanglement generation. We 
also use the behaviour of (F) to determine if the accumulative noise is a product of 
local noise of the type introduced in section [5] 

7.3.1. Decay noise (low-T limit) - Earlier in section r7.2.1l we found that using the EB 
as a single-qubit channel with local decay noise results in the accumulative noise 
being independent of N. Thus all chain lengths were equivalent to a chain with just one 
spin. When the EB is used to mediate a c-cr^ gate with local decay noise, our numerical 
results show that (F), for the lengths investigated, is independent on N. However, 
further investigation reveals that A^^l itself does possess a weak length dependence 
demonstrating that this result is a consequence of {F) being completely insensitive 
to these changes. A length dependence of Al^l is expected since, in contrast to the 
single-qubit channel in section 17.2.11 the projected master equation in the subspace 
Ti-o © Til ffi does not have commuting coherent and dissipative contributions. Our 
numerical calculations furthermore indicate that Emin is weakly length dependent but 
is never positive over the parameter range and lengths investigated. 

We find that the coupling at which the fidelity drops to (F) < 0.99 is Kf/J = 
4.0 X 10""^ independent of N. In addition to the average gate fidelity we also 
compute the specific gate fidelity of the initial state | ++) using equation ([3]) 

and the numerically determined superoperator A^^l. Using this initial state makes the 
operation equivalent to the noisy generation of a two-qubit cluster state. We find that 
the coupling at which this fidelity drops to F++ < 0.99 is k'^^ / J = 8.5 x lO^'^, and 
turns out to be independent of N. Thus, this specific preparation is twice as resilient 
to decay noise than the average preparation. 

7.3.2. Dephasing noise - For local dephasing noise we find that emin, shown in 
figure llOf a). is never positive over the parameter regime considered and therefore 
the ability to generate entanglement is retained in the presence of this noise. In a 
similar way to the single-qubit channel fidelity we find that the gate fidelity, plotted 
in figure [TOT b). rapidly converges with increasing N. To understand the nature of 
the accumulative noise we attempted to fit the fidelity to a model where the ideal 
two-qubit gate U is implemented and then product noise Aj^'^^ (g) A|^|^^ is applied 
where Aj^j^^ is a single-spin superoperator describing noise from the class introduced 

in section [S] In fact we found that the best fit was obtained when A|^|^j was further 
restricted to the case used in section [7.2.21 where 7 = k and /3 = (k. The validity of 
these fits determined from Fa had a peak infidelity of 1 — Fa ~ 0.3 for strong coupling 
with N = 2 and so the actual noise bears no resemblance to this product model. This 
rapidly drops to 1 — Fa ~ 5 x 10^^ for larger N, which indicates that the product 
noise model becomes more valid for longer chains. This behavior is sensible since in 
long chains the two spin packets do not overlap for the majority of the evolution time 
T and therefore experience independent noise during this time. We also find that the 
coupling at which the fidelities (F) and F++ drop to 99% are k//J = 2.5 x 10^^ and 




Figure 10. For local dephasing noise - (a) The minimum eigenvalue Emin of the 
partial transposition of the 4 qubit mixed state for the bipartition (12) (34) 
and (b) the average gate fidelity {F), both as a function of the chain length A'^ 
and coupling strength k/J. 



k'j^ /J = 6.5 X 10 ^ respectively, for chain lengths up to TV = 12. 

7.3.3. Thermal noise (finite-T) - For local thermal noise we find that emin, shown 
in figure [TTT a). is only marginally increased and remains negative over the parameters 
we considered. Despite this from our analysis of the single-qubit channel there is 
good reason to suspect that for longer chains and higher temperatures this noise 
will generate Emin > 0. Along with (F) depicted in figure [TlTb') Emm has a very 
weak dependence on N for T/A < 0.2 similar to that encountered for a single-qubit 
channel in section 17.2.31 We again fit the fidelity surface with a product noise model 
^mod assuming the same single spin noise as used in section 17.2.31 where a{T) and 
/3(r) remain unchanged from those in equation ([8]) and ^{T) — Ca{T) with a fit 
parameter This fitting produces a linear dependence of with N as found earlier 
for the single qubit channel. In this case the infidelity for the product noise model 
is 1 — Fa < 2.3 x 10~^ and is therefore a good approximation to the accumulative 
noise over the parameter range investigated. This indicates the noise is very effective 
at eliminating correlations which might be built up by the dynamics of the chain. We 
also find that the temperatures at which the fidelities (F) and F++ drop to 99% of 
their T = value are T//A = 0.29 and T^+/A = 0.34 respectively, for chain lengths 
up to TV = 12. 

7.3.4- Depolarizing noise (high-T limit) - For local depolarizing noise we observe in 
figure I12f a) that emin becomes positive for a sizable portion of the parameter range 
explored. We can therefore only be certain that entanglement generation is possible 
with the EB outside this region of parameters, i.e. k/J < 0.15 for chain lengths up to 
= 12. The average gate fidelity plotted in figure [T2r b) shows that k ^ J is required 
in order to achieve a reasonable average fidelity. When fitting a product noise model 
to {F) we obtain an infidelity 1 — Fa < 4.5 x 10~^ which decreases significantly with 
larger coupling k/J. This is consistent with stronger local noise destroying correlations 
between the spin packets and therefore decorrelating the noise. We also find that the 
coupling at which the fidelities (F) and F"*""*" drop to 99% are Kf/J = 8.0 x 10~^ and 
nt^ /J = 1.4 X 10^'^ respectively, for chain lengths up to A = 12. 
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Figure 11. For local thermal reservoirs - (a) The minimum eigenvalue eniin of 
the partial transposition of the 4 qubit mixed state for the bipartition (12) (34) 
and (b) the average gate fidelity (_F), both as a function of the chain length A'^ 
and temperature fcgT/A using a coupling k/J = 0.02 




Figure 12. For local depolarizing noise - (a) The minimum eigenvalue Eniin of 
the partial transposition of the 4 qubit mixed state for the bipartition (12) (34) 
and (b) the average gate fidelity (F), both as a function of the chain length A'^ 
and coupling strength k/J. The critical coupling Kc/</ for each chain length A'^ 
at which the entanglement generating capability of the EB is no longer certain is 
given by the intersection with the £^\^ = plane shown in (a). 



8. Implications for graph state generation 

In section [3] a general scheme for creating arbitrary graph states was outlined whicli 
exploits the multi-qubit circuit implemented by the EB. A complete characterization 
of the influence of noise on the full scheme is beyond the scope of the current work. 
Here, we instead focus on the most direct implication of the results presented in section 
17.31 by investigating graph state generation for five qubits using the EB to mediate a 
two-qubit c-cr^ only. As an example we focus on the generation of a linear cluster state, 
shown in figure [TST a). and a GHZ state, shown in figure ITST b). under the influence 
of local dephasing noise. These two states were chosen because they have the same 
number of edges but very different topologies. 

The generation protocol begins by initializing the first five qubits of a ten qubit 
register in the state | +) . For the previous calculations in section 17.31 the input and 
output qubits were exclusively the end spins in the EB. For graph state generation, 
however, each usage of the EB to establish an edge will necessarily involve using 
different spins of the chain as the input qubits (as well as output qubits). Our 
numerical calculations indicate that the superoperator A^^' does depend weakly on 
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Figure 13. We consider the generation of (a) tlie linear cluster state and (b) 
the GHZ state for five qubits. The qubits representing the graph vertices are 
initialized in the first half of the register as shown in (c). Two schemes for the 
sequence in which the EB is used to mediate the necessary two-qubit gates are 
investigated. Scheme (i) proceeds with the gates in a sequential order, while 
scheme (ii) proceeds by performing gates between qubits which are closest in the 
register at a given step. For the linear cluster state the two schemes are explicitly 
shown in (c). In (d) the graph state fidelity Fg for the linear cluster state ('o') 
and GHZ state ('□') constructed with scheme (i) is plotted, and in (e) the same 
is shown for scheme (ii). 



the input qubit locations. For this reason the order in which the gates are performed 
will affect the quality of the overall state. To illustrate effect this we computed the 
graph state generation protocol for two different gate ordering schemes. The first, 
scheme («), performs the necessary gates in a simple sequential order according to the 
qubit labels. The second, scheme (ii), performs a gate between qubits, from the list 
of edges to be established, which are closest in the register at a given step (which 
is often not unique). These two schemes are shown explicitly for the linear cluster 
state in figure fTSFc). We computed numerically the superoperator A^^' for all input 
locations required for these two schemes to generate the linear cluster state and the 
GHZ state. By concatenating the appropriate noisy gates in the specified order the 
imperfect graph state was obtained. 

For scheme (j), shown in figure [T^f d). we find that the linear cluster state fidelity 
Fg is slightly above that for the GHZ state with them having a fidelity greater than 
99% for k/J< 5.7 x 10^* and k/J < 6.0 x 10^"*, respectively. For scheme (ii), shown 
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in figure [T5f e) , the GHZ state fidelity is virtually unchanged from that of scheme (i) , 
whereas the linear cluster state fidelity has dropped below that of the GHZ state giving 
a slightly reduced n/J < 6.1 x 10~^ for a fidelity greater than 99% to be attained. 
The average gate fidelity for a iV = 10 chain, computed in section [7751 using the end 
spins as the input qubits, drops to 99% when k/J < 2.6 x 10"'^ for local dephasing 
noise. Thus the noise tolerance for the generation of these graph states is reduced 
by a factor of approximately 4.5 in comparison. Additionally, the results appear to 
indicate that only a weak change in the fidelity is observed with the graph topology 
and input qubit locations for a TV = 10 chain. Further work is need to confirm if this 
inscnsitivity is maintained for more varied topologies over larger numbers of vertices 
and longer chains. 

These preliminary results give a clear indication of useful directions for future 
work. This includes studying the dependence of the EB on the input qubit locations 
with different noise models as well as chain lengths and determining wether there 
is any generic behaviour. For the full graph state generation scheme of section [3] 
the quality of the multi-qubit circuits implemented by the EB with noise needs to be 
studied along with their likely dependence on the input qubit locations. This may well 
reveal a trade-off between using the EB less times with more qubits or more times 
with fewer qubits. With this information and for a given graph topology, the graph 
state generation scheme under the influence of noise could be optimized, in terms of 
the sequence and type of gates implemented. 

9. Conclusion 

In summary we have investigated the influence of local noise on the mirror inverting 
EB spin chain. For the case where the EB is utilized as a single qubit channel we have 
found that the accumulative noise of the EB for local decay noise is independent of 
its length N and explained this unexpected behaviour. Additionally, we have found 
that neither local decay nor dephasing noise cause the EB to become entanglement 
breaking, in contrast to local thermal and depolarizing noise. For the latter two 
cases we have determined the critical length and critical coupling /tc/ J at which 
entanglement breaking occurs, respectively. The local depolarizing noise Kc/ J is found 
to exhibit a power-law dependence on TV which is explained by the competition between 
the speed and spreading of the spin packet in the EB. 

For the case where the EB is used to mediate an entangling c-cr^ gate we flnd 
that the entanglement generating capability of the EB is never lost in presence of 
local decay or dephasing noise. For local dephasing noise the resulting operation of 
the EB becomes progressively closer to product noise with increasing N . Both local 
thermal and depolarizing noise are well approximated by product noise models due to 
the severity at which they decohere the spin chain. For local depolarizing noise our 
results indicate that the entanglement generating capability of the EB can only be 
guaranteed for couplings k/ J < 0.15 for chain lengths up to = 12. 

We have also performed a preliminary analysis of the graph state generation 
scheme which uses the EB, with local dephasing noise, to individually implement the 
necessary c-a^ gates for a five qubit graph state. As an example we focussed on the 
linear cluster state and GHZ state. As expected we found that the concatenation of 
these operations reduces the tolerance to noise. In this case the reduction was by a 
factor of approximately 4.5 for both states compared to the average fidelity of a single 
gate. More work is needed to determine how this behaviour scales for larger graph 
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states constructed with longer EB chains. Finally, we note that while our results 
indicate that there are tight constraints on the levels of tolerable noise for graph state 
generation this can be weakened within the framework of one-way quantum computing 
due to the separation between the preparation and consumption of entanglement |39j . 
In principle this enables the resulting noisy graph states to be purified prior to their 
use gOlIll]. 
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Appendix A. Limits of the local noise model 

In section [5] the general form for the Kraus operators was given in equation ([7]) for the 
noise model considered here. Since the expressions for the general matrix elements 
of the diagonal operators are rather lengthy we have postponed their introduction to 
this appendix. To begin with the following quantities are defined [7] 

AoW = J[l + 2e-5("+/3+27)* + e-("+«*], 
\i{t) = \2{t) = \[l-e-(»+f^)% 
Xs{t) = i[l - 2e-5("+/3+27)* + 

using the definitions from section [5l We then find that for /i 7^ 



Ti = ^{x - 2/.t - Ao + Aa)^ 



T2 - i(2/i + 2;-Ao + A3)^ 



T3 = i(2/i + x + Ao- Ag)^ 



T4 = i(a; - 2^ + Ao - Ag)^ 





Ao 


+ X3 - X 




4/x2- 


(Ao- 


A3) (a: - Ao 


+ A3) 




Ao 


+ A3 - a; 




4/x2- 


(Ao- 


A3) (a; - Ao 


+ A3) 




X -i 


- Ao + A3 




4^2 + 


(Ao- 


A3) (a; + Ao 


-A3) 




X -i 


- Ao + A3 




4/^2 + 


(Ao- 


A3) (a; + Ao 


-A3) 



making one further definition x — ^ 4/i2 + (Aq — A3)2. As expected the Kraus 
operators in equation ([7|) simplify considerably for a number of important special 
cases. In particular for /i = the Kraus operators become E\ — -v/Aq 1, E2 = -\/A]"ct^, 
— \/A2 fT^ , £"4 = a/As ct^ . For (3 = k and a = 7 = they reduce to 

^'-^ Vp-i ' ' " I 



with pi — e describing pure decay noise; for 7 = k and a — (3 — Q they reduce 
to El — y^l, E2 — VI — P2 with p2 — \J 5(1 + e"***) describing pure dephasing; 
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while ior a — [3 ~ "f — they become Ei — ^1 — 3p3 1, E2 — ^/pscr^ , E3 — y/pacr^, 
E4 — ^Jpia^ with p3 = -1(1 — e""*) describing depolarizing noise [3T]. For finite-T 
noise with 7 — and a + /5 = k the Kraus operators take the form 



\/e^ \ ^ rrr ( 



For /i 7^ the average fidelity for a single-qubit channel experiencing this class of 
noise can be computed from equation Q as 

1 /^o , (Ao-A3 + a;)2(Ao + A3 + a;) , (Aq - A3 - a;)^|Ao + A; 



3 



6 V V + (Ao - A3)(Ao - A3 + IV - (Ao - A3) (A3 - Ao + 
For finite-T noise with 7 = this equation reduces to 
1 1 1 



2 3 6 

and has a double exponential decay independent of T making it equally applicable to 
the T = decay channel. For the special case where ^ = the expression for (F) 
dramatically simplifies to 

(^^,,=0) = \ (4Ao + 2) , 

from which the well known fidelities |31| for dephasing and depolarizing channels can 
easily be evaluated as 

(-Pdophasc) ^3^3^ ' 



1 1 



(-Fdcpolar) — „ + r. 



-Kt 



2 2 

which exponentially decay to | and i respectively. 



Appendix B. Positive partial transposition criterion 

Detecting the presence of entanglement in a general bipartite mixed state can be 
achieved to an extent through the use of the positive partial transposition (PPT) 
criterion [34l [35] . For an operator O acting on the Hilbert space H, = C^" ® C^*" of 
two systems a and h the partial transpose with respect to system a is defined as 

da db 

0'^" = ^ ^ (i,m|0| j,7i) I j,TO) (i,7i| , 

— \ mn—1 

in terms of some basis {| z, m) | i — 1, • • • , da, m = 1, ■ ■ ■ : df,} of H. Although this 
definition is basis-dependent the spectrum of O^"^ is not. Note also that while the 
transposition of the full system ab preserves the positivity of the full density matrix 
{^pTayr^ — -> general the transposition with respect to any subsystem does 

not, and so the partial transpose is not a completely positive operation. 

A state p of the system ab is separable if and only if it can be expressed as a 
convex combination of product states 

X 
i=l 
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with Pi > and J2f=iPi — 1- The PPT criterion then states that p^"^ > is a 
necessary condition for separability of any da x db systems. Importantly for 2 x 2 and 
2x3 systems the PPT criterion is necessary and sufficient for separability [331 133 ■ 

Appendix C. Jamiolkowski isomorphism 

In section |4] extensive use of the Jamiolkowski isomorphism |361 [7] is made. This 
isomorphism establishes an equivalence between quantum states and superoperators. 
To begin suppose we have a type of subsystem S with Hilbert space Ti. of dimension 
d and spanned by basis states {\i) | i = 0, • • • , d — 1}. For such a subsystem any 
density operator p can be expanded in the operator basis {| i) (j | | i, j = 0, • • • , d — 1} 
with its corresponding matrix elements contained in a d^-dimensional vector pij. 
Superoperators are defined as linear, trace-preserving, completely positive maps of 
density operators to density operators. Consequently, a superoperator A^-'^l acting on 
a subsystem 5* is completely described by a d^ x d'^ super-matrix with elements in the 
operator basis A^^'j.^ as 

aW{N)0-|}- EAW,jfc)(n. 

This information can be mapped to a quantum state by using two copies a and 
b of the subsystem S initially prepared in the maximally entangled state | ^I*) = 
J2f=o I ® ^^"^ applying the superoperator A'^l to b as 

(1 ® aW){| VP) 1} = ^ xi I 0' I ® 0' 1} = 

ij=0 

as shown in figure [3]^a) . The resulting density operator p^ for the two S subsystems 
then completely describes A^-'^l by containing all of its operator matrix elements 
A[^1{| i) {j |}. This can then be used to compute Al^ljp'"} = on any state via 

A[ii{p'"} - d xi ( E f^^i Pkijj I (j I - /^°"*' 

ij=0 \ kl=0 ) 

thereby giving the inverse isomorphism. 

For a superoperator A^^l acting on two S subsystems the Jamiolkowski 
isomorphism proceeds in an analogous way by applying A^^l on one half of two 
maximally entangled pairs of subsystems with the setup depicted in figure EJb) as 

(la ® t^ll ® (* ® I *) (VI/ ^ pA. 

The resulting density matrix for the four S subsystems is then 
^^-^2 E l^>(Jl®APl{|^fc)(J7|}®|fc)(^|, 

ijfei=0 

and again completely describes A^^l through its matrix elements in the operator basis 
of two subsystems. 
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